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Abstract 

We extend an earlier, configuration space method to find the Wilson coefficients of operators ap- 
pearing in the short distance expansion of thermal correlation functions of different quark bilinears. 
Considering all the different correlation functions, there arise, up to dimension four, two new op- 
erators, in addition to the two appearing already in the vacuum correlation functions. They would 
contribute substantially to the QCD sum rules, when the temperature is not too low. 



The operator product expansion of two currents at short distance |f[ played a key role in the 
past in writing down the extensively used QCD sum rules The vacuum matrix elements of the 
local operators appearing in the expansion were used to parametrize the nonperturbative effects 
of the QCD theory at low energy. Using dispersion representations for the correlation functions, 
these sum rules relate the multitude of low energy physical quantities, like the masses and coupling 
constants of hadrons, to these parameters of theoretical interest. 

In extending these sum rules to finite temperature, two sources of complications arise. One 
is the interaction of the current with the particles of the medium (heat bath), modifying the 
dispersion representation Q. The other is the breakdown of Lorentz invariance by the choice of 
the reference frame, in which the heat bath is at rest at a definite temperature. The residual 
0(3) invariance allows, in general, more operators of the same dimension in the operator product 
expansion than at zero temperature 

While the effect of the medium was included in the very first paper § on thermal QCD sum 
rules, the non-perturbative contributions due to the additional operators were seriously considered 
much later || . So far only the thermal sum rule for the vector current has been written with these 
contributions and analysed at low temperature It still remains to explore the consequences of 
various sum rules for correlation functions of different quark bilinears. One ingredient necessary to 
write them is the set of coefficients with which the operators, at least the leading ones, contribute 
to the different sum rules. 

Earlier we employed a configuration space method to calculate the coefficients of gluon op- 
erators up to dimension six, which appear in the vacuum expectation value of correlation functions 
of quark bilinears §,§. It is based on finding the singular solution for the quark propagator in 
the background field of the gluons. It finds a natural extension to thermal correlation functions. 
Here we find the coefficients of all operators up to dimension four, appearing in all correlation 
functions built out of any two quark bilinears. We start with a brief review of this method. 

Let i/j(z) be the quark field of a single flavour of mass m, where its SU(3) colour indices are 
suppressed. Its equation of motion is 

(-ip+m)^{z) = 0, D li = d ll -iA li , A il = gA i il \ i l% (I) 

where A fJi (z) is the colour gauge field and A 1 are the SU(3) Gell-Mann matrices. Assume, to begin 
with, that A^ is a given classical field. Then the time ordered product of the quark fields may be 
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split as, 

Tip a {z)i) b {y) = -iS a b(z,y)+ : ^ a {z)4>b{y) :, (2) 

where a, b are the spinor indices. The first term is a c-number function, containing the short 
distance singularities of the time ordered product and satisfies 

(-ilp z + m)S(z,y) = 5\z-y). (3) 

The second term is a bilocal operator with finite matrix elements as z — > y and involves quark 
dynamics at large distances. 

We now outline the method of solving eq.(3). We set^m = 0. It is very convenient to choose the 
Fock-Schwinger gauge fL0| |, (z — y) fi A^(z) = 0, where the gauge potential A^(z) can be expanded 
in terms of the field strength and its covariant derivatives at z = y. In the following we put y = 
and get, 

A^z) = \z a G ail {Q) + ■ ■ ■ . 

Then the solution will automatically be gauge covariant. We convert (3) into a second order 
differential equation by setting 

S(z,0) = ipE(z,0). (4) 

Then E(z, 0) satisfies 

(a-P^-Q)E(z,Q) = 5 4 (z), (5) 

where 

P tl (z) = iz a G atl + ---, Q{z) = pGj+---. 

Here and below we use matrix notation, eg., 7G7 = 7^ l G ' iiv r f '. The fundamental or the singular 
solution to (5) is now given by the Hadamard Ansatz for i?(z,0)^| O], 

E(z,0) = ^1 + V(z)ln(4z 2 fi 2 ) + W{z)j , (6) 

where V(z) and W(z) are analytic functions of z, regular in the neighbourhood of z = 0. Although 
not singular, W(z) is added to ensure that E(z,0), by itself, satisfies eq.(5). The factor 4/^ 2 in the 
argument of the logarithm represents a mass scale (not related to the quark mass), arbitrary at 
this stage. Inserting this Ansatz in (5) we get the equations for V and W, which can be readily 
solved in power series in Zp. 

The solution for S(z,0) is worked out in ref.(8) up to gauge field strengths of dimension six. 
Here we reproduce this solution up to dimension four, 



S(z,0) 



__L zGG7 + _L (zG7)(7G7 ) + ^{\lGG 1 - -i-G- 1-} j /»N:V 
i^GG,- 1 -^?} 



(7) 



Again the last term, which is regular, ensures that S(z,0), by itself, satisfies (3) up to dimension 
four. Here we have omitted terms involving the gauge field covariant derivatives, as they do not 
survive the trace over colour in the correlation functions. In this case S(z,0) satisfies, 5(0, z) = 
S(-z,0) 

Once we have found the c-number coefficients multiplying the classical gauge field strengths 
in S(z,0), the latter can be considered as gauge field operators of QCD theory. Eqs. (2) and (7) 



2 Operators linear in m can be easily sorted out later. 

3 For the time ordered propagator z 2 stands for z 2 — it 
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thus constitute the short distance expansion of the two point function of the quark fields. Clearly, 
higher order corrections to the coefficients due to gluon propagation is beyond the present method. 

Consider now the short distance expansion of the time ordered product of two quark bilinear 
operators, 

Ji(s) =: ${z)T X il>{z) :, J 2 (z) =: ^(z)T^(z) :, 

where I\ and T 2 are any two of the Dirac matrices. Using (2) it can be written as a sum of three 
types of terms, 



Tj 1 (0)j 2 (o) = 4f ) +4t ) +4 re9) 



"12 ' 



where 



t%p = trT 1 S(z,0)T 2 S(0,z), (8) 
fit } = -«(r 2 S ro (0, z^U : M0)Mz) ■ -ipiSmiz, 0)T 2 ) ab : Mz)M®) :, (9) 
t<" ff) = :#)r^(^(0)r^(0):, (10) 

where tr indicates trace over both Dirac and colour matrices. In (9) we have included the missing 
term in S(z, 0), which is linear in the quark mass, 

7 777, 

S m (z,0) = S(z,0) 



4ir 2 z 2 ' 



The first term gives the coefficients of the gluon operators. The second term, after expanding ip(z) 
in the bilocal operators as, 

1>(z) = V(0) + *"D^(0) + • • • , 
gives the coefficients of all quark operators, 



t 



1 



{i(r 2 ^i - r^r 2 ) ab + mz 2 (rir 2 + r 2 ri) a6 /2} : v> a V>& ■. 



w = 

12 2tT 2 (z 2 ) 2 

+(r 2 *r 1 + r^Uz" ■. ^D^ b ■] (11) 

The third term is of no interest to us. 

Thus 4? an d t^) constitute the most general operator product expansion of two quark bi- 
linears, giving the coefficients of all the gluon and quark operators up to dimension four, having 
various tensor and spinor structures. Clearly at the level of operators there is no reference to 
temperature. It is in taking the matrix elements of these operators that the difference between 
the vacuum state and the thermal state shows up. In taking the vacuum expectation value, the 
different tensor and the spinor operators are projected on to their corresponding Lorentz scalars. 
The situation is different for the thermal expectation value, which for any operator O is defined as 

< O > T = Tre-P H 0/Tre-P H , 

where H is the QCD Hamiltonian, (3 is the inverse temperature T and the Tr(ace) is over any 
complete set of states. Here the choice of the thermal rest frame, where the temperature is defined, 
breaks the Lorentz invariance and we have to include 0(3) invariant operators as well. 

It is, however, convenient to restore Lorentz invariance formally in the thermal field theory by 
introducing the four-vector the velocity of the heat bath. [ u 2 = 1 and = (1,0,0,0) in the 
rest frame of the heat bath.] Thus while at zero temperature scalar operators are formed from 
tensors by contracting its indices among themselves, at finite temperature additional scalars can 
be formed by contracting its indices with u M . 

Consider first the piece t^p . The Lorentz covariance at finite temperature now allows us to 
write 

< tr c G a pG\ a > T = {gaXgpa ~ 9aa9f)\)A 

- {u a u x g fia - u a u„gp\ - upu\g aa + u fj u a g a \)B. (12) 
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By contracting indices on both sides, we get 

A = ^<G a ap G a ^>T+l<u a Qi f ^> T , (13) 
B = \<u a Q a ap u > T , (14) 

®«/3 being the traceless, gluonic part of the stress-tensor of the QCD theory, 

Q 9 af3 = ~G a aX G^ + \g aP Gl a G^. (15) 

Since we are interested in the coefficients of operators of dimension four only, we may already use 
(12) in the expression for S(z,0), which then simplifies to 

°) = i U^) i^rh'G^ + |(o* + 14)) , (16) 



— (u ■ zf - 2/n(4zV), b = 8u ■ zln{Az 2 ^ 2 ) 



where 

4 

(G) 

Insert this expression for S(z,0) in t^' and use (12) again for the dimension four operators 
generated in the product. Then it only remains to evaluate the traces over 7-matrices for different 
correlation functions (specified by T\ and T 2 ) to find the coefficients of the operators in A and B. 

Next we turn to the piece t^). As with the gluon operators, we project the quark operators 
on to the corresponding scalars. At finite temperature the most general decomposition of the two 
operators in (11) in spinor space, conserving parity, are given by 

<Mb>T = ^<^>T+^<#V>T (17) 
< TpaiD^b > T = (j^baU + U^) ba V + U^SbaY + dj^baZ (18) 



where 



u = ^ < w > T ~ < u"e£„u" > T , v = \< «"e£X > T , (19) 



j _ jjj _ | _ jjj _ 

Y=-<^iu-D^> T +—<^> T , Z= — <^iu-D^> T -—<^ii)> T (20) 
12 12 12 

and 0£„ is the traceless, fermionic part of the stress tensor, 

_ l 

Now ipiu ■ Dip can be reduced to mip^ip and a total derivative. Further < tpi/iip >t= 0, for zero 
chemical potential and we shall not consider it any further. We are thus left with only < ipip >t 
in (17) and U and V in (18). Inserting these decompositions in and evaluating the traces over 
the 7-matrices, we can get the coefficients of these operators. 

We have thus reduced the calculation of the Wilson coefficients of operators appearing in the 
expansion of the product of any two quark bilinears to evaluating traces over 7-matrices. These 
coefficients have to be Fourier transformed to momentum space for use in the QCD sum rules. 
Here we only state this result for the correlation function of vector currents (Fi = 7^, T 2 = 7„). It 
can be expressed through two invariant amplitudes, T\ and T 2 , say, which arc functions of uu = u ■ q 
and q 2 , 



if. 



:<"'-■ <r./ /( (-.)./,,(<V) r {-r 1 ^ + ^jT 1 (cj,q 2 )+^-^J \ u v - ^ J T 2 (u; , g a ) 

(21) 
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Then we find (Q 2 = -q 2 ), 



2m n 2 R / 111! 2 " 

t x = - ¥ <^> T+ ^ ¥ <G^> T - w (i-^)<uef u > T 

g 2 \ ( 2w 2 \ , 2 2 . 5 7w 2 
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{ (l - ^) ( ln (Q 2 /^ 2 ) + 27) - | + ^} < n9% > Tl (22) 

T 2 = < „e'« > T (WW) + 2 7 -£) <«e 9 o T , (23) 

where 7 is Euler constant. 

The large q 2 behaviour of the Wilson coefficients can be improved by the renormalisation group 
equation, ft effectively shifts the renormalisation scale of these coefficients from /i to y/Q 2 . The 
operators g 2 G il _ v G^ v and mtptp and hence their coefficients are invariant under this change of scale. 
But it mixes the operators 0^„ and 0^„, a situation familiar from the analysis of deep inelasic 
scattering Let us write the contributions of Q$ and 9 to either Ti or T 2 , for short, as 



C f n f Q% + C g G 



y 1 



where the coefficients Cf >g [= C f, g (log(Q 2 / /j, 2 ) , g([i 2 ) , to 2 / q 2 )] can be read off from (22, 23) and n/ 
is the number of flavours. The coefficients satisfy the coupled renormalisation group equations 
with an anomalous dimension matrix |l3| . Diagonalising this matrix, we see that the combination 
of the coefficients (rifCf + ^C g ), corresponding to the total energy- momentum tensor operator 
(rifQf + 9 ), has zero anomalous dimension, while the combination [Ct — C g ), corresponding to 

the operator n/(^0^ — 9 ) has anomalous dimension, — /f w ^ d, with d = f + nf). Thus 

C f n f Qf\, + C g Q% 

1 f, ^ 16=, w „ f =, w 16 



(n f C f + yC s )(n/0^ + 0») + an,(C, - C s )(y9^ - , (24) 



where C /l9 = C/, a (l,p 2 (Q 2 )) and a = {(g 2 (^ 2 )/g 2 (Q 2 )}- d/2b , b = 11 - 2n//3. Thus ignoring 
g 2 -corrections (which for C/ is not calculable anyway by the present method), the stress tensor 
contributes to Ti, say, asQ, 

Ti - - (16 + 3 n/)g2 (l - 2^ 2 /9 2 ) {> + o(16e//3 - e s )| M } , (25) 

where we use notation for the rest frame of the heat bath, e = Uf€f +e g , with ej and e g being the 
energy densities of quarks and gluons respectively. 

Finally let us list the set of thermally averaged operators contributing to the thermal QCD 
sum rules. Considering all the different correlation functions and restricting to operators up to 
dimension four, there are two new operators, namely, e = (rifCf + e g ) and (16e//3 — e g ). Two 
condensates in the vacuum sum rules, now appearing as thermal averages, viz, < ijjtp >t and 
< G'^ l/ G^ a >t, have their temperature dependence given by jbij, 

dP 

< ipip > T =< 0\iptp\0 > — — , 

am 

< ^G%G^ a > T =< 0\£^G%G^ a \0 > ~(e - BP), 

where P is the pressure. In the hadronic phase they are known from chiral perturbation theory 
Jl5| , |l6| ] and at high enough temperature from ordinary QCD perturbation theory. Thus only one 

4 Our result does not agree with the one stated by the authors of the second paper in ref. (6), because they work 
in the deep inelastic limit and do not use the renormalisation group equation. 
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condensate, viz, (16e//3 — e g ) cannot be readily evaluated. It may be obtained from one of the sum 
rules themselves. It should be noted that though the contributions of the new operators would be 
small for small enough temperature, they would, in general, dominate over the condensates in the 
QCD sum rules in the neighbourhood of the phase transition, where the latter ones begin to 'melt' 
away. 

I have much pleasure to thank H. Lcutwyler for many helpful discussions. I also thank B.L. 
Ioffe for bringing the references (6) to my attention. I gratefully acknowledge the hospitality at 
the University of Bern, where this work was carried out. 
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